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Abstract 

We analyse the evolution of a quantum oscillator in a finite temperature environ- 
ment using the quantum state diffusion (QSD) picture. Following a treatment similar 
to that of reference we identify stationary solutions of the corresponding Ito equa- 
tion. We prove their global stability and compute typical time scales characterizing 
the localization process. The recovery of the density matrix in approximately diag- 
onal form enables us to verify the approach to thermal equilibrium in the long time 
limit and we comment on the connection between QSD and the decoherent histories 
approach. 



1 Introduction 

One of the approaches to quantum theory, developed the last years, is the quantum state 
diffusion (QSD) picture [||, ||, |^. The construction of this picture was motivated by stochas- 
tic reduction theories and the necessity of describing individual experimental outcomes. 
Mathematically it is equivalent to the Lindblad master equation (hence it applies only to 
the Markovian regime of open quantum systems). Its main feature is that it describes the 
stochastic evolution of an individual system in Hilbert space, a treatment complementary 
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to the deterministic evolution provided by the density matrix. Being a phenomenological 
picture, it has proved within its domain of apphcabihty to be in very good agreement with 
experiments involving individual systems |Q and it is helping us to obtain more physical 
intuition about processes of the microworld. 

Whenever the master equation for the reduced dencity matrix p describing the evolution 
of the open system is of Lindblad form there is a unique and consistent way of unravelling 
it into an ensemble of individual system state vectors each of which obeys a stochastic 
differential equation. Thus if the evolution of p is given by, 

f = ~[H, P] + Y^i^L^pLi - hlKp - ipLt L„) (1) 

the evolution for the state vector of the individual system is given by the non-linear 
stochastic Ito differential equation, 

\d^) = '-H\i^)dt + Y.{{Li).^K-\LiL^-U^^^ (2) 

Here H is the Hamiltonian of the system in the absence of environment, sometimes mod- 
ified by terms depending on the Lindblad operators L„, which model the effect of the 
environment. The independent stochastic differentials d^n are complex and satisfy, 

Mdin = 0, M{dindU) = 0, M{dindCJ = 5mndt (3) 

In the above ()^ denotes mean value with respect to and M expresses the average over 
the ensemble. The density matrix is recovered by the formula 

p = M(|V^)(^|) (4) 

and consistency between the two pictures means that when mean values taken with respect 
to \ip) in QSD are averaged over the ensemble, results agree with those obtained using p. 

2 The Model 

In this paper we apply the QSD picture to the class of models consisting of an harmonic 
oscillator in a finite temperature dissipative environment. The Hamiltonian of the system 
in the absence of the environment is given by, 

H = huj{a^a + ^) (5) 

while time evolution is described by a Lindblad master equation (1) with two environment 
operators, 

Li = [(n + l)7]i/2«, L2 = {n^f'^a^ (6) 
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Here a and are the annihilation and creation operators for the harmonic oscillator, 

a = {(Tp/h)q + i{aq/h)p (7) 

where cxg and Up denote the dispersions of position and momentum for a standard coherent 
state, 

and 

n=[e-x.Y>{hu/kT)-l\-^ (9) 

while the coefficient 7 represents dissipation. The process under study could be the damp- 
ing of a mode of the electromagnetic field in a cavity, coupled to a beam of two-level atoms, 
some of which are initially excited or the damping of a coherent light beam propagating 
in a weakly absorbing medium [^]. 

The treatment will be similar to that of reference 0. There, the QSD picture was 
applied to quantum Brownian motion (QBM) models resulting in approximate stationary 
solutions to Ito equation. These were correlated coherent states. They minimize a more 
general uncertainty relation and are localized in phase space. For quadratic potentials 
the solutions are exact and stability is global i.e. every initial state evolving under Ito 
takes the form of the stationary one after some time. One of the main results of this work 
was that, after localization time, every initial density matrix tends to the form, 

P = J fiP,q,'t)\lPpg){lPpg\ (10) 

The above expression gives a natural way for approximately diagonalizing p in the same 
basis at all moments of time, expressing phase space localization in the density matrix 
language. It resembles the well known P-representation of the density matrix ^ but differs 
in that \ippq) are correlated rather than standard coherent states and f{p,q,t) is always 
positive by construction. The expected approach to thermalequilibrium was also proved 
in Ref. [0] and finally localization rates were studied and the connection with decoherent 
histories was exemplified using phase space projectors constructed from \ippg). 

The linearity of Lindblad operators in position q and momentum p makes our analysis 
resemble the one of QBM model but the physics is obviously not the same. Here we are 
studying a process, Markovian at any finite temperature T (even small) and not only in 
the Fokker-Plank limit (which is a high temperature limit), studied in ref. [0. The effect 
of the environment here is incorporated into two Lindblad operators expressing different 
processes and the sources of phase space localization together with the localization rates 
turn out to be of a different nature in both cases. 



The same system has also been subject of study previously [|T^, [Tl|, [T^. In ||12| the 



aproach to thermal equilibrium was shown numerically, while in [10| it was argued that we 



have coherent states as stationary solutions to Ito equation and a localization theorem was 



proved in [|rT|. Our study concentrates on diffferent aspects of the problem (and we also 



believe that the authors of W^ have made a mistake in their calculation). The purpose 
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of this paper is to study connections between localization and decoherence as in ref. 
but in a model different to that one, and also to give an analytic account of the numerical 
results of ref. [0. 



3 Stationary Solutions 



Following [7] in the search for the stationary solution we require that they satisfy the 
condition, 

X % % 

+ |#) = exp[-gc/p- ^pdq + -d(j)]\ijj) (11) 



with. 



^tdt + Y.{cPrMn + <P:dC) (12) 



n=l 
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dq = Ldt - (7/2)g + ^ [a{q, Ln)d^n + a(L„, q)dQ (14) 

Condition (10) means that the shape of the stationary solution is preserved and only p and 
q change. Here is a phase and by p and q we mean (p)^ and (q)^ respectivelly. We have 
also introduced: 

a(f,6) = (ft6)-(ft)(6) (15) 
After lengthy but not difficult calculations condition (9) leads to, 

{q\p q) = ij{q) = exp[-A(g - qf + -pq] (16) 

with A real, 

which is a standard (harmonic oscillator) coherent state characterized by minimum uncer- 
tainty: 

cry, = h'/A (18) 
This property makes the stationary solutions look like phase space points to classical eyes. 



.2 - 



qdt - (7/2)p+ ^[(T(p,L„)rfe« + a(L„,p)dC] (13) 



n=l 

2 



4 Localization Rates 

The existence of stationary solutions gives rise to the questions of stability and localization 
rates. To prove that every state tends to the form of the stationary one, we observe that 
the stationary solution is eigenstate of a. To show global stability, we need to show that, 

Md{^af < (19) 
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Now for a non Hermitian operator O we can define the spread as, 

(A6)2 = (t(6, 6) (20) 
Using the two above equations we obtain, 

2 2 

Md{Aaf = ^{n + l/2)[h^ - AR^ - 2^(Ap)^ - 2^{Aqf] (21) 
where R, is the symmetrized correlation between q and p, 

R^l{{p,q})-{p){q) (22) 

while (Ap)^ and (Ag)^ denote the dispersions. Then setting, 

{Aqr = al{l + Q) {Apr = alii + P) (23) 

clearly 

P > -1, Q > -1 (24) 

and stationary solution is obtained for P = Q = R = 0. The dispersion of a reads in terms 
of P and Q, 

(Ad)' = (P + Q)/4>0 (25) 
which substituted in (20) leads to 

MdiAaf = -27(n + 1/2) + ^ + ^ + ( (26) 

The above expression is negative and vanishes only for a coherent state. This proves the 
global stability of the stationary solutions. Our expression for the rate of localization is 
not in agreement with the one obtained in ||TT|. 



It is obvious from (25) that the localization process is characterized by a timescale of 
order ti = [7(?T+ l/2)]~^. This coresponds to a minimum rate of localization, 

^rf(A«)2 ^ ^ ^/2)(Ad)2 (27) 

Then using expression (8) we obtain, 

~ - tanh(huj/2kT) (28) 
7 

which for huj << kT gives. 
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while for huj » kT we have, 

t/ ~ - (30) 

7 

When studying the approach to thermal equilibrium the meaning of (30) will become 
evident. In the case of an initial state consisting of a superposition of wavepackets a large 
distance d apart, (Ag)^ ~ (P and term will be the dominant one in (25). Then in the 
high temperature limit, we obtain by virtue of (22) and (24), 

t ~ (^-^^ 
d'^rwykT 



This result agrees with the usual "decoherence time" for high T [16 



5 Thermal Equilibrium 

The approach to equilibrium in Ref. was shown by computing numericaly the behaviour 



of (n) and the occupancy probabilities for various number states \n). Using our stationary 
solution we may do the same analitically.For times greater than the localization time, 
{n(t))^ is given by, 

2 2 

(n(t))^ ^ {nlt))st = ^f + ^p' = a*{t)a{t) (32) 

Here a denotes the eigenvalue of a. Using equations (13) and (14) we may derive an 
evolution equation for a{t) |T3[| , 

a{t) = a(0) exp[-{iu + j/2)t) + y/Wy f' exp[-{iuj + 7/2)(t - t')]d^2{t') (33) 

Jo 

It is clear that in the long time limit only the part of the solution with the integral over 
the stochastic process will survive. Then (n(t)) reads, 

{n{t)) = n^ fexp[{^ + iu;)T + i^-tiuy]dUr)dC2ir'). (34) 
JO 2 / 

From the above expression only when we take the mean over the ensemble we are able to 
say that once equilibrium is reached (n) will fluctuate around n as expected. Then the 
standard result, 

M{h{t))st = n (35) 

is recovered, since M[d^2{'r)dC,2{T')] = 6{t — r')dTdT'. Using our stationary solution, in 
principle, we are able to compute time averages over the stochastic process and thus repro- 



duce the numerical result of Ref. [|T2| for the occupation probabilities analytically. (This 
result is needed to exhibit the thermal nature of the fluctuations in (34)). However, the 
manipulations such a computation would involve are cumbersome. We shall then proceed 
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using the density operator to treat this problem. We may recover p following the treatment 
of Ref 0. Once localization has taken place Eq° (4) is of the form, 

p = j fip,q,'t)\p(i){pq\ (36) 

Hence every initial p approaches the phase space diagonal form (36) on the localization 
timescale. Equation (36) is the P-representation of p with the nice property that f{p,q,t) 
is positive. A Fokker-Plank equation is easily derived for f{p,q,t) |[T3[| : 

—f = -^—f + - — (qf) + mu^q—f + - — (pf) + I^^f + !^n^muj—f (37) 
dt mdq 2 dq dp 2 dp 2mu d(f 2 dp^ 

It admits a unique Gaussian stationary solution being globally stable and approached for 
times t >> 

/s(p,g) = 7r^exp(-=|ap) (38) 
iT^n n 

This is the P-symbol expression for a thermal density operator. Hence every initial state 
approaches thermal p for long times, as expected. Then the elements {n\p\n) may be 
evaluated and give the standard result, 

(nlpk) = =(l + rr^^)-("+^) (39) 
n 

for the occupation probabilities, since the probability of observing n quanta in a coherent 
state(2.e. |(^|m)P) is given by the poissonian distribution with mean |ap. We note that 
our result together with that of Ref. (which is numerical and for a specific initial state) 
give the first traces of a proof that the ergodic hypothesis holds for our system. 

One of the featurs of our model is that it is applicable at low temperatures as well. When 
the temperature T ^ 0, we see from equations (8), (33) and (36) that the system decays 
to its ground state as expected, (for zero temperature fs{p,q) becomes infinitely sharp 
around a = as seen from (38). A very interesting result obtained in connection with 
equation (30) is that localization and relaxation to thermal equilibrium proceed essentialy 
on the same time scale for T —>■ 0. This should be expected since for n = 0, Li = 'ja 
and 1/2 = 0, expressing only dissipation of energy to the environment. Therefore in this 
limit the mechanism of decoherence is not very efficient, exhibiting the deeper connection 
between environment aly induced decoherence and noise. 



6 Comments 



In the context of explaining the emergence of classical behaviour as a result of the interac- 
tion of a system with its environment, the approaches of "decoherence of density operators" 



|T^, |T3|, |T7| and of "decoherent histories" [|T^, [2^, ^ give the frameworks for a 
systematic treatment. In the former the approximate diagonilization of the reduced den- 
sity operator, expressed by (36) for this class of models is essential while in the latter the 
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decoherence of histories is the prerequisite to emergent classicahty. The close connection 



between QSD and decoherent histories pSf can be illustrated using the results obtained so 
far. The localization properties of the stationary solutions to Ito equation makes natural 
the study of phase space histories. We thus take phase space projectors of the form, 

Ppg = / dpdq \pq){pq\ (40) 

Here Tpg is a phase space cell. Obviously those projectors are only approximate ones, the 
validity of approximation depending on the area covered by the cell and the nature of its 
boundary |2^. Construct now histories with the above projectors and assume our initial 



state is p, then the decoherence functional is given by, 

D{a,g!_) = rr(P,„K,t J^a„_,....<n^ai<HPo]Pa'J....P<_jPaO (41) 

Here a, 0/ denote the two strings of projections at succesive moments of time 
characterizes the phase space cell at a time moment and Kl^__^ is the propagator of 
the density operator obeying the master equation of our model. D{a,gf_) given by (42) is 
recovered from the decoherence functional of the closed system after having traced out the 
environment and under the assumptions that the initial density operator factorizes and 
the process is Markovian. 

For intervals between projections longer than the localization time the density operator 
will always evolve to the diagonal form (37). Studying then (42) we conclude that because 
of the diagonality of p the off- diagonal elements D{a,gf_) will be approximatelly equal to 
zero. Therefore approximate decoherence is achieved. Then probabilities computed are 
found to be most strongly peaked around the classical path and are essentialy the same 
for both approaches provided we always refer to scales grater than |. This is because the 
probabilities of phase space histories (computed from the diagonal elements of D) and the 
probabilities in QSD are both equivalent to the probabilities derived from the description of 
our system using the Fokker-Planck equation (38). Thus all the conclusions of the analysis 
of Ref. [0 are valid here as well. However, a comment we would like to make is that in 
our model localization and decoherence occur in the Lindblad operator while in Ref. ^ do 
not. This is due to the fact that here the eigenstates of L are preserved by the Hamiltonian 
while in general this need not be the case. 

In connection with the results obtained in we can then conclude that the diagonality 
of p in a basis consisting of eigenstates of the same operator at all momements of time, 
the construction of decohering phase space histories (requirments for emergent classicality) 
and the study of the approach to thermal equilibrium may be acomplished in a natural way 
using the stationary solutions to Ito equation. Hence the description of processes in terms 
of individual systems, the localization properties of its solutions and the close connection 
with decoherent histories, make quantum state diffusion very useful both on intuitive and 
calculational grounds in the study of the emergence of classical behaviour of open quantum 
systems. 
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